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The neutral current effects of the future high statistics atmospheric neutrino 
data can be used to distinguish the mechanisms between a oscillation to a 
tau neutrino or to a sterile neutrino. However, if neutrinos possess large di- 
agonal and/or transition magnetic moments, the neutrino magnetic moments 
can contribute to the neutral current effects which can be studied by the single 
7T° production events in the Super-K data. This effect should be included in 
the future analyses of atmospheric data in the determination of to tau or 
sterile neutrino oscillation. 
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I. INTRODUCTION 



Neutrinos might possess two properties which are feeble, but will be important barom- 
eters of physics beyond the standard model scale. These are neutrino mixings (masses and 
oscillations) and neutrino magnetic moments JJ. Before 1998, experiments gave bounds 
on these properties, in general. 

But the recent results from the Super-Kamiokande Collaboration |3[] have provided a 
strong evidence for a deficit in the flux of atmospheric neutrinos, which are presented in the 
form of the double ratio 

(NJN e ) MC ' k ' 

which implies the existence of oscillation. The measured value of R for Super-Kamiokande 
is 0.61 ± 0.06 ± 0.05 for the sub-GeV data and 0.67 ± 0.06 ± 0.08 for the multi-GeV data, 
while we expect R = 1 in a world without oscillations. Muon neutrino oscillation into another 
species of neutrino provides a natural explanation for the deficit and even the zenith angle 
dependence. The — + v T oscillation is the most favorable solution for the atmospheric 
neutrino problem, whereas the — > v e oscillation is strongly disfavored by CHOOZ results 
||. The oscillations into sterile neutrinos [y s ) give a plausible solution as well ||. This 
evidence for the neutrino oscillations is also supported by the SOUNDAN2 || and by the 
Super-Kamiokande [pj and MACRO || data on upward-going muons. 

It is usually assumed that the neutral current effect in neutrino oscillation experiment is 
unchanged, since any standard model neutrino produced by oscillation has the same neutral 
current (NC) interaction. Therefore, the ratio of the neutral (NC) and charged currents 
(CC) events is important to investigate the neutrino neutral current. Thus the observation 
of single ir° events, induced by the neutral current, by Super-Kamiokande Collaboration || 
can lead to important physical implications. 

Since the n° NC event is detected as two diffuse rings whereas the CC events due to v e 
are detected as one diffuse ring due to e ± and one sharp ring due to n ± and the CC events 
due to Vp are detected as two sharp rings from ^ and 7r ± ||. Thus, a NC event can be 
discriminated from a v e CC event and a CC event [^,|TT|]. It has been considered difficult 
to separate NC and CC events clearly. But the single 7r° events described above can be 
used to discriminate the NC events from CC events. Indeed, it is believed that the cleanest 
way to identify NC events in Super-Kamiokande is to detect a single 7r° from the process 
v + N — > v + N + tt°, with iV being either a neutron or a proton below Cerenkov threshold. 
The 7T° is detected via its decay into two photons which lead to two diffuse e-like rings whose 
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invariant mass is consistent with the ir° mass ||. The ratio of 7r°-like events to e-like events 
compared to the same ratio of the Monte Carlo in the absence of the oscillation has been 
measured by the Super-Kamiokande Collaboration 0, 

Rir°/e = (7r ° /e)data = 0.93 ± 0.07 stat ± 0.19 sys , (2) 
(ti - u /e)MC 

where the systematic error is dominated by the poorly known single tt° cross section, and 
the statistical error is based on 535 days of running. The ratio R n o/ e is expected to be 1 for 
— v T oscillations while 0.75 for — v s oscillations or — v e oscillations if one takes the 
measured z/ M /z/ e ratio to be 0.65 [|Kj. The admixture of v^, — v T and — v e oscillations also 
leads to a deviation of R n o/ e from 1. Therefore, a precise measurement of the ratio will be 
used to distinguish — > v T from — > v s oscillation. 

At first sight, it is likely that any deviation of R n o/ e from 1 implies muon neutrino 
oscillation into a sterile neutrino. However, if there exists a large muon neutrino magnetic 
moment (diagonal or transition), it will produce an additional neutral current effect which 
has to be separated out to draw a definite conclusion. Indeed, right after the discovery of 
the neutral current, the upper bound on the muon neutrino magnetic moment was given 
EJ. Also, the experimental bounds on transition magnetic moments and other properties 



in view of NC data were presented [|T3 



The theoretical problem of obtaining a large neutrino magnetic moment has begun with 
interactions beyond the standard model 0. In general, the loop diagram will have a (mass) 2 
suppression, presumably by M x where Mx can be the W boson mass or a scalar mass. It 
is possible to have a large Dirac neutrino magnetic moment if the loop contains a heavy 
fermion 0, 

(3) 



Ml 



where m is the mass of the heavy fermion. This mechanism can be generalized in models 
with scalars fTS . 



But, the same loop without the external photon line would give a contribution to the 
neutrino mass matrix. Therefore, one expects, taking the coupling as 10~ 3 , 



^~ 10 'w x {mf" B ~ n (4) 

where \x& = eh/2m e c is the electron Bohr magneton and we used m u ~ 0(1) eV for the 
numerical illustration. To suppress the contribution to the mass and still allow a large 
magnetic moment, continuous |L5| and discrete |l6j symmetries have been considered. In 
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this case, the neutrino magnetic moment can be as large as \i v ~ 10~ 10 fiB, which is not 
affected by the SN1987A constraint ~ lCr 13 /!^ |L7j since this bound applies to the electron 
neutrino only. 

However, a large transition magnetic moment to a sterile neutrino is not forbidden that 
severely. For example, one can introduce a transition moment with an accompanying mass 
as large as several hundred MeV. Of course, the masses of the light neutrinos are bounded 
by eV. In this case, the transition neutrino magnetic moments can be as large as 10~ 7 /xg and 
may contribute to NC events. In particular, we are interested in the single 7r° production 
through a large transition magnetic moment, which would contribute to R^o/e- In this spirit, 
we will obtain the lower bound on the transition neutrino magnetic moment (to a sterile 
neutrino) from R^o/ e . 

This paper is organized as follows: In Sec. 2, we describe the amplitude for the single 
7T° production. In Sec. 3, the kinematics and the differential cross section for the single ir° 
production are given. In section 4, we present the contribution to the cross section of the 
7T° production generated by a possible neutrino transition magnetic moment. In Sec. 5, we 
discuss the physical implications based on the numerical result. 

II. PRODUCTION OF THE SINGLE NEUTRAL PION 

The single ir° production has two contributions: one from the production and decay of 
the (3/2,3/2) baryon resonances and the other from the continuum contribution. At low 
energies (E < 2 GeV), the contribution from baryon resonance production is a dominant 
one for a single 7r° production ]Tj| : 



v + N -> u + N*, 

N* -> vr° + N, (5) 

where N* represents baryon resonances. The cross section for the single ir° production in the 
region W < 1.6 GeV (W is the hadronic invariant mass in the final state) can be described 
following Fogli and Nardulli |]18||. The effective Lagrangian for the neutrino neutral current 



is defined by 

Ave = ^^(1 + 75Wf °, (6) 
by assuming the following general V, A structure of the hadronic part of the NC: 

^ C = 9 3 vV? + g\A\ + g*V* + g\A\ + gyV® + g\Al (7) 
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where A\{i = 3, 8, 0) are the SU(3) nonet partners of the CC Neglecting the strange 
and charm NC's, one can write 

J?° = 9vV x 3 + g A A\ + g' v V? + g' A A®, (8) 

where 

^° = y|w + ^), (9) 

A'^^Al + V^Al), (10) 
the electromagnetic current being given by Jf 71 = + l^' . From the Weinberg-Salam 



model [19 



gv = ^-sm 2 9 w , 9A=~, (11) 
^ = -~sin 2 w , g' A = 0. (12) 



The isospin decomposition of one 7r° channels is given by 

A{yp ^ vw *) = i(2A 3 + A x ) + ^S, (13) 

A{yn ^ vn ^) = I(2A 3 + A x ) - J\s. (14) 
The reduced matrix elements A\, A% are given by 

A3=^(A° A + Al + A N ), (15) 

Ax = ^A° NNn - V2AI - ^=A° N + A° s + A° P + A° D , (16) 

where A®, % = N, NNn, P, S, D are given in the appendix and the indices S, P, D denote 
Su, Pn, Du | ]2U[| , respectively. The contributions to the amplitude S come from the following 
amplitudes 

S=\V3 U\a% + ^ + \a% + \ A % + \A^ . (17) 

As is well known, the dominant contribution to the amplitude A comes from the A reso- 
nance in this region (W < 1.6 GeV) [^TJ. Then, the amplitude for the single ir° production 
can be described by 

A NC ~ ^l a J Q (18) 
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where 



l a = u(k') la (l + j 5 )u(k), (19) 



-(g^@-Q^) l5 g v £L 



■(g^Q ■ p' - Q^p'^gy^ - g^g A C£ 



M N 

v 

Al 



u(p), (20) 



where Mn is the nucleoli mass, D w is the propagator of Rarita-Schwinger field which is 
given by 

j) + M' / 2 Pfl p p l pyn P -Pulp 1 \ rgn 
u wW) - p2 _ M , 2 + . MT [9w 3 M , 2 + 3 M/ 3 7m7 P; | , ^-U 

and q = p' + q n — P = k — k' is the momentum transfer, T is the decay width, Cj and 
Cf(i = 3, 4, 5) are the vector and axial vector transition form factors as defined by Llewellyn- 
Smith p2]| , and M' is mass of the A resonance. As shown in Ref. [|11J, the form factors 
Ci(q 2 ) can be obtained by comparison with the values of the helicity amplitudes given by 



the relativistic quark model [E3J. The explicit forms are given by 



1.7y/l - q 2 /AM 2 R 



p i„i\ y 1 i n (22) 

3Kq J [1 -gV(M fi + M^) 2 ] 3 / 2 [l-g 2 /0.71 GeV 2 ] 2 ' 1 ; 

C^) = -^C^q% (23) 

C,(q 2 ) = 0. (24) 

Then, the vector form factors CY{q 2 ) used in Eq.(20) are given by 

CY{q 2 ) = V3Q(q 2 ). (25) 

The axial form factors Cf(q 2 ) is also taken to be the general formula 

CV) = irSw (26) 

where M A = 0.65 GeV. 

III. KINEMATICS 

Consider the process given in (5). It is convenient to choose the center of momentum 
frame. Without loss of generality, we can choose the initial four momenta of the neutrino 
and the nucleon in the CM frame as (p,p, 0, 0) and (En, —p, 0, 0), respectively, where 
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V 
E N 



1 + 2E u /m N 
m N + E v 



(27) 
(28) 



'1 + 2E u /m N 

where E v is the incident neutrino energy in the LAB frame. The final four momenta of 
neutrino, nucleon and 7r° are k',p' and q n , respectively, where 



k' = E v , [1, (cos 6, sin 9,0)] , 
p' = E N . 



(29) 



1, 



\ 



m 



— ^L, 

AT' 



— cos (3 ■ cos 9 + sin /3 • sin 9 ■ cos 0, 



— cos (3 ■ sin 9 — sin /3 • cos • cos 0, — sin (3 ■ sin 0) 
<2V = -Ett 



(30) 



L 



mi 



, 1 — -y(— cos a • cos — sin a • sin 6* • cos0, 



cos ck • sin 9 + sin a ■ cos ■ cos <f>, sin a • sin <f>) 



where 



cos a 



cos (3 = 



2\k'\-\q v \ 
k' 2 +jf 2 — q w 



(31) 

(32) 
(33) 



2\k'\ ■ \f\ 

The angles 9, <p correspond to the rotations around z— and x— axes, respectively Then, the 
differential cross section da can be expressed as 



da = 



(2n) A \A 



NC 



4(p • k) 



JVC 



32(27r) 4 M JV S I/ 



dE u /dE n d(pd(cos 9) , 



where 



with 



Q2 

\Anc\ — —^-L^J^, 



51 fyvi 



spins 



and 



J = - V J 



(34) 



(35) 



(36) 



(37) 



spins 

where the summation is performed over the hadronic spins. 
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IV. CONTRIBUTION FROM THE NEUTRINO MAGNETIC MOMENT 



In this section, we consider the A production arising from the Feynman diagram shown 
in Fig. 1. The decay of A to a nucleon plus 7r° follows in the detector, and one observes the 
v' + N + 7T° final state. The — v 1 — 7 vertex is parametrized by a transition magnetic 
moment 

ifv'v^Buil^uXJ^U^)^ (38) 

where q = I — I' = p' + q n — p is the momentum transfer. The coupling f v i v at q 2 = is 
the transition neutrino magnetic moment in units of the electron Bohr magneton and will 
be denoted as /'. 

The squared matrix element that describes the single 7r° production, induced by a tran- 
sition neutrino magnetic moment /', can be written as 

\A M \ 2 = f^-M^J™, (39) 

where 

M^ = \ £ [u(l^ a q a u(l)]mK^Mn] (40) 
spins 

and 

= \ E J e r j r (4i) 

spins 

where 3 e ™ is the hadronic electromagnetic current given before by J^ m = + (1/3)V^°. 

V. RESULTS AND DISCUSSIONS 

In this section we present the numerical results of the cross section of the single n° 
production for the NC interactions. The calculated cross section generated by the neutrino 
transition magnetic moment is shown in Fig. 2 as a function of the incident neutrino energy 
for the hadronic invariant mass less than 1.6 GeV. We find that the values of the cross section 
is of the order KT 40 cm 2 for E v < 2 GeV. 

In order to see how the contribution to the cross section generated by the neutrino 
magnetic moment can be constrained by the experimental results of the ratio R n o/ e , it is 
sufficient to calculate the ratio -^wu where af,(jQ C are the cross sections of the single 7r° 
production from the neutrino magnetic moment and the standard model NC interactions, 
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respectively. The reason is that e-like events (CC) is not affected by the presence of the 
neutrino magnetic moment. In Fig. 3, we plot this ratio 



as a function of the incident neutrino energy E u for /' = f' Q = 0.6 x 1CT 8 . f' G is defined 
as the value giving a similar contribution as the NC interaction. If /' is e times f' Q , Fig. 3 
should be multiplied by a factor e 2 . Note that the contributions from the transition magnetic 
moment and from the standard model NC do not mix in the process is + N — > v' + N' due the 
unmixable 7 matrix structure among these two. However, for u + N — > u' + A there are terms 
which mix these two contributions. This is because the form factors T^, defined in AT^N, 
can match the 7 matrices through taking out one index by q^. The dashed line corresponds 
to the case of no cut whereas the solid line corresponds to the hadronic invariant mass cut at 
1.6 GeV. The current experimental result of the ratio R n o/ e implies that the possible excess 
from 1 amounts to 0.13 from which we can obtain the constraint on the neutrino magnetic 
moment. For example, at E v = 5 GeV a constraint r/v/wc < 0.13 leads to 



The transition magnetic moment of this magnitude implies the muon neutrino and sterile 
neutrino mass matrix of the form 



where mi 2 ~ m 2 i is roughly 10 4 times mu ~ O(10~ 2 ) eV. Thus mi 2 is of order < 100 eV. 
The diagonalization process should not change the mass of the muon neutrino drastically, 
i.e. > ^i2/ m n ~ 1 MeV. Therefore, a singlet neutrino at the intermediate scale with 
possible beyond the standard model interactions (scalar or gauge) can lead to a sizable 
transition magnetic moment. 

The effects of the muon neutrino transition magnetic moment should be separated out 
toward a final determination of the muon neutrino oscillation to the tau neutrino or to a 
sterile neutrino. The most promising method is to study the energy distribution of the final 
7T°, since the kinematics for the magnetic moment is different from the NC interactions where 
the former has a 1/g 2 dependence in the differential cross section while the latter has no q 2 
dependence at low energy. 

Indeed, if the transition magnetic moment is discovered by the measurement of the energy 
distribution, it will hint an intermediate scale physics. On the other hand, one can compare 



<7 fi 

r f'/NC = -jvc 



(42) 



/' < 2.2 x 10 



-9 



(43) 




(44) 
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this anticipation to the earlier expectation that neutrinos must oscillate due to the belief 
that singlet fermions, the remnants of grand unification or the standard model superstring, 



would be present at the intermediate scale p4j , f2q| . Similarly, if singlet neutrinos are present 
much above the eV scale, there may be a large transition magnetic moment which can be 
detected by the future high statistics atmospheric neutrino experiments. 
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APPENDIX 



We present the expressions of A®, i = N, NNtt, P, S, D |18|| : 



A% = -?>^-G F J x V2g NN ^u( V ') 



x 



F 2 , . 

g'v F il\ + 9v7}Jj-[ 7A, Q\ ~ 9Ao F Al\l5 



N 



1 

3' 



p" + $+M N 



A° NN « = ^g f j x V2 9nn Mp')t^ 



N 



k 2 -M N 



F 2 1 



u(p), 



A P 



■-J-G F J x f P u(p') ^ _ m2 D 6 A g' Alxl5 u(p), 



A° s = 9 -^J X fsu(p')^ 
s 2 V 2 k — M 



Ml 

g' v G s xlx - <&-^g-[ 7 A, ^]7s + 3a^7a 



'2Mat 



«(p), 



4° 



(45) 



(46) 



(47) 
(48) 



(49) 



($ S pA - QVkyT^ + (fc ■ Q 9 pX ~ Q p k x )g{ 

M N 



H 



V M 2 



N 



(p " Q ^ - " 9 pX 9'a\h s a ^ 



u(p) 



(50) 



where G F = 1.023 • 10~ 5 M^ 2 is the Fermi constant, Qnn-k/^ = 14.8 is the NNir coupling 
constant, Mm is the nucleon mass, m„ is the pion mass and Mr the mass of the generic 
resonance R. We also assume Ma = 0.65 GeV/c 2 (Ma = 1-0 GeV/c 2 for the axial mass 
of the P 33 (Pn, Sii, Dis) . The generic vector form factor in the amplitudes Eqs.(45)-(47) is 
given by the general formula 



V(t) = V p (t) + V n (t), 



(51) 



where V p (V n ) is the electromagnetic form factor with a proton (neutron) as target. 
With regard to the pion and nucleon form factors, we use the following usual forms, 

1 



F?(t) 



1 - t/0.47 GeV 2 ' 
3.7t \ 



1 



1 



4M 



N 



t) (l-t/0.71 GeV 2 ) 2 ' 



(52) 
(53) 
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F 2 y (t) = 1.855 1 T - ; 5-, (54) 

2W V AM n) (1 -t/0.71 GeV 2 ) 2 ' V ' 

F A (t) = 1.23- :— (55) 

w (1 - 1/0.81 GeV ) 2 v ; 

where t = Q 2 is the moment transfer and the axial form factor is characterized by Fa(0) = 
1.23, which is derived from neutron decay, whereas Ma = 0.90 GeV/ c 2 is perfectly compatible 
with an overall fit to neutrino experiments. The form factors G{ and ifj are explicitly given 
in Ref. |B|. 

The axial form factors in the amplitudes Eqs.(45)-(47) are taken to be the general formula 



A(0) 
'l-t/Mlf 



Mt) = , t (56) 



where Ma = 0.65 GeV for and M A = 1.0 GeV for the other resonances. 
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FIGURES 




FIG. 1. Feyman diagram for the A production arising from the neutrino transition magnetic 
moment. 



15 








4 



1 



E v (GeV) 

FIG. 2. The cross section generated by the neutrino transition magnetic moment as a function 
of the incident neutrino energy for the hadronic invariant mass less than 1.6 GeV. /' is taken to be 
0.6 x 1(T 8 . 
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f = 0.6 x 10" 



Without Cut 
W< 1.6 GeV 
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FIG. 3. Plots of Tfi /nc as a function of the incident neutrino energy for /' = 0.6 x 10 8 . The 
dashed line corresponds to the case of no cut and the solid line corresponds to the invariant mass 
cut at 1.6 GeV 
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